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Abstract 


DUNIC,  RONALD  L.  A  Method  for  Correcting  and  Reducing  the  Pressure 
Data  Recorded  at  the  PAM  II  Stations  During  GALE.  (Under  the  direction 
of  Allen  J.  Riordan.) 

^  The  data  recorded  at  the  PAM  II  stations  during  GALE  will  enable 
researchers  to  gain  a  better  understanding  of  the  mesoscale  processes 
involved  in  the  various  phenomena  ,t-hat  occurred-  dwitvg-OAL^r  however, 
the  recorded  pressure  cannot  be  directly  analyzed  because  it  has  not 
been  corrected  or  reduced  to  sea  level.  This  paper  presents  a  method 

for  correcting  the  pressure  data>  and  then  reducing  the  corrected 

6.  ,  - 

pressures  to  sea  level. 

.  _,The  Barnes  objective  analysis  scheme  was  employed  to  determine  the 
pressure  correction  for  each  PAM  II  station.  The  corrected  pressures 

are  reduced  to  sea  level  using  the  standard  NWS  reduction  procedure 

..  •  .*  ‘ '  Wc*  S 

outlined  in  the  Manual  of  Barometry  (1963),  This  method  proved  to  be 
-•rather  Involved  and  include^a  number  of  assumptions  to  obtain  the  sea- 
level  pressures.  An  empirical  formula  was  derived  for  the  GALE  dataset 
to  calculate  the  sea-level  pressures  for  the  PAM  II  stations.  Pressure 
analyses  involving  the  corrected  and  reduced  PAM  II  pressures  for  two 
cases  were  in  good  agreement  with  the  NMC  analyses  for  the  same  times. 

Another  method  for  representing  the  surface  geostrophic  wind  field 
and  the  horizontal  pressure  fiel<t  introduced  by  Pielke  and  Cram 
(1987),  is  explored  to  determine  the  true  ageostrophic  wind  components 
in  a  cold-air  damming  situation  and  to  produce  a  horizontal  pressure  — 


field  that  would  eliminate  the  arbitrariness  of  the  standard  pressure 
reduction  procedure.  The  result  is  an  improvement  over  the  geostrophlc 
wind,  but  the  observed  winds  are  still  highly  ageostrophic  In  this 
cold-air  damming  situation.  '  -  f 
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1.  INTRODUCTION 


Weather  forecasters  continue  to  be  frustrated  by  unforecasted 
weather  events.  Every  forecaster  has  experienced  the  surprise 
thunderstorm  or  unexpected  intensification  of  a  storm  system  causing 
public  scrutiny  of  their  forecasting  ability.  In  some  cases  the  result 
may  be  nothing  more  than  a  picnic  being  rained  out,  but  in  others, 
property  damage  and  even  loss  of  life  may  result. 

Table  1  shows  the  spatial  and  temporal  ranges  of  the  synoptic  and 
meso-  scales.  Synoptic-scale  organization  has  been  the  basis  for  much 
of  the  success  in  current  weather  prediction;  however,  many  damaging  or 
even  life-threatening  phenomena,  such  as  heatwaves  and  droughts,  arctic 
outbreaks  and  heavy  snows,  flash  floods,  tornadoes,  hail, 
thunderstorms,  windstorms  of  hurricane  force,  and  downbursts  from 
thunderstorms  are  all  mesoscale  phenomena.  Many  times  these  events  are 
not  predicted  because  the  forecaster  is  unable  to  resolve  the  mesoscale 
processes  involved  due  to  the  inability  of  the  current  National  Weather 
Service  (NWS)  surface  and  upper-air  station  network  to  adequately 
resolve  weather  systems  below  the  synoptic  scale.  Thus,  the  current 
weather  charts  do  not  show  everything  the  forecaster  needs  to  see  in 
order  to  provide  a  better  prediction. 

Another  difficulty  in  forecasting  is  that  the  mesoscale  processes 
involved  are  not  well  understood.  What  is  needed  is  a  more  dense  data 
collection  network  that  would  enable  forecasters  and  researchers  to 
identify  and  study  mesoscale  features  as  they  develop  and  propagate  to 
gain  a  better  understanding  of  the  processes  at  work  and  subsequently 
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Table  1.  Spatial  and 

temporal  ranges  for 

the  synoptic  and  meso-  scales 

(Dirks,  et  al. ,  1988) . 

Scale 

Spatial 

Temporal 

Synoptic 

>  2000  km 

24-48  hrs 

Meso-a 

200-2000  km 

3-24  hrs 

Meso-p 

20-200  km 

<  3  hrs 

Meso-T 

2-20  km 

<  3  hrs 

produce  a  better  forecast. 

In  1973,  the  development  of  a  computerized  measurement  system  to 
be  used  in  the  study  of  a  wide  range  of  mesoscale  phenomena  including 
hail  and  thunderstorms,  squall  lines,  sea  breezes,  and  regional  air 
pollution  was  begun.  This  system,  called  Portable  Automated  Mesonet 
(PAM),  collects  and  displays  mesoscale  meteorological  data  from  a 
surface  array  of  remote  sampling  systems  (Brock  and  Govind,  1977).  The 
data  recorded  at  the  PAM  stations  include  temperature,  pressure, 
humidity,  wind  speed  and  direction,  and  rainfall  amounts  (Brock  and 
Govind,  1977).  This  data,  when  incorporated  into  the  National  Weather 
Service  (NWS)  station  network,  enables  researchers  to  study  the 
processes  involved  in  mesoscale  storm  systems. 

The  first  operational  use  of  a  PAM  network  was  during  the  National 
Hail  Research  Experiment  (NHRE)  conducted  in  northeast  Colorado  in 
1976.  During  this  experiment,  the  PAM  network  demonstrated  its  ability 
to  "now-cast"  surface  mesoscale  patterns  associated  with  thunderstorm 
development.  The  experiment  also  showed  that  the  PAM  network  offered 
researchers  a  reliable  source  of  continuous  mesoscale  data  as  routine 
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service  on  the  PAM  stations  was  not  required  (Brock  and  Govind,  1977). 

The  Genesis  of  Atlantic  Lows  Experiment  (GALE)  conducted 
during  the  winter  of  1986  deployed  a  network  of  fifty  PAM  II  stations 
through  North  and  South  Carolina  and  Southern  Virginia  in  order  to 
gather  data  to  enable  researchers  to  study  the  mesoscale  processes 
associated  with  East  Coast  storms.  These  PAM  II  stations  measured 
in-situ  temperature,  wet-bulb  temperature,  station  pressure,  wind  speed 
and  direction,  humidity,  dew  point  temperature,  and  precipitation 
amounts. 

Of  the  meteorological  data  recorded  at  these  and  all  weather 
stations  around  the  world,  the  only  parameter  that  is  not  directly 
analyzed  is  the  station  pressure.  If  one  constructed  a  chart  with 
isobars  of  station  pressure,  the  result  would  look  more  like  a 
topographic  map  rather  than  a  useful  weather  map  because  pressure  is 
strongly  dependent  on  elevation.  In  order  to  create  more  useful 
weather  charts,  the  station  pressure  is  often  reduced  to  a  standard 
reference  level  to  eliminate  the  effects  of  elevation  on  the  pressure 
field  and  allow  forecasters  to  see  the  true  pressure  pattern. 

By  convention,  the  reference  level  chosen  is  sea  level  because  it 
is  assumed  that  pressure  reduced  to  sea  level  will  enable  forecasters 
to  determine  the  surface  circulation  pattern  and  horizontal  pressure- 
gradient  forces  (Sangster,  1987).  Unfortunately,  finding  a 
satisfactory  method  for  reducing  the  station  pressures  to  sea  level  has 
not  proven  to  be  an  easy  task  for  meteorologists  over  the  past  century. 
The  difficulty  has  involved  finding  a  representative  way  of  depicting 
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the  temperature  of  the  fictitious  column  of  air  between  the  station 
elevation  and  sea  level. 

Due  to  the  concentrated  network  used  during  GALE,  the  pressure 
data  recorded  at  the  PAM  II  stations  may  permit  a  true  mesoscale 
definition  of  the  pressure  gradient.  From  the  pressure-gradient  field 
the  horizontal  pressure -gradient  force  and  the  geostrophic  wind  can  be 
uniquely  determined.  However,  before  the  PAM  II  data  can  be  reduced  to 
sea  level  the  data  must  first  be  corrected  to  enable  the  PAM  II  network 
to  fit  into  the  NWS  station  network. 

The  research  described  in  this  thesis  deals  primarily  with 
enabling  the  pressure  data  recorded  at  the  PAM  II  stations  to  be 
incorporated  into  mesoscale  analyses  and  studying  a  method  of  better 
representing  a  surface  pressure  field  from  which  surface  geostrophic 
winds  can  be  computed. 

1.1  Literature  Review 

In  order  for  the  PAM  II  data  to  be  utilized  in  resolving  the 
mesoscale  features,  it  is  important  for  the  data  to  be  meteorologically 
real  and  not  artificially  Induced  by  sensor  or  measurement  errors 
(Wade,  1987).  One  method  for  checking  data  quality  is  to  use  reference 
sensors  transported  to  each  site  as  was  done  during  GALE.  The  method 
provides  a  useful  check  on  the  accuracy  of  the  sensors  as  long  as  the 
calibration  of  the  reference  sensors  is  maintained  and  a  sufficient 
number  of  checks  are  accomplished  so  as  to  provide  a  statistically 
significant  correction  (Pike,  1985).  Due  to  the  remoteness  of  some  of 
the  PAM  II  stations,  there  was  not  a  sufficient  number  of  checks  during 
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the  GALE  Project. 

The  errors  commonly  observed  in  mesonet  data  can  result  from 
failure  of  a  component  in  the  system,  calibration  errors,  sensor  drift, 
exposure  of  sensor  to  the  elements,  or  system  noise  (Wade,  1987).  The 
detection  of  these  errors  requires  documentation  of  the  sensor's 
performance  compared  with  other  sensors  in  the  network.  This  can  be 
accomplished  through  a  direct  one-to-one  comparison  with  other  stations 
or  through  an  objective  analysis  (Wade,  1987).  In  this  research,  buddy 
checks,  comparisons  between  two  nearby  stations,  were  performed  with  a 
handful  of  PAM  II  stations  and  their  nearest  NWS  station  to  determine 
the  feasibility  of  employing  an  objective  analysis  to  obtain  the 
corrections.  As  was  done  in  the  quality  control  approach  of  Wade 
(1987)  this  research  employs  the  Barnes  objective  analysis  scheme  to 
obtain  the  corrections  to  the  PAM  II  pressure  data. 

In  addition  to  sensor  error,  the  offsets  obtained  from  the  Barnes 
analysis  and  buddy  check  methods  are  influenced  by  real  sub-grid  scale 
events  such  as  thunderstorms.  For  the  most  part  these  sub-grid  scale 
events  are  short-lived  and  the  variability,  or  noise,  caused  by  these 
events  is  smoothed  by  averaging  over  a  number  of  cases. 

Once  the  pressure  data  have  been  corrected,  a  method  of  reducing 
these  corrected  pressures  to  a  common  reference  level,  usually  sea 
level,  can  be  found.  As  was  stated  earlier,  this  has  not  proven  to  be 
an  easy  task.  During  the  first  thirty  years  of  its  existence,  the 
Weather  Service  used  no  fewer  than  six  different  methods  for  the 
reduction  of  pressure  to  sea  level,  with  each  method  improving  upon  its 
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predecessor  until  the  Weather  Service  was  confident  it  had  found  a 
useful  algorithm.  A  historical  review  of  the  early  methods  and  the 
problems  associated  with  them  is  provided  in  Appendix  8.1. 

The  Weather  Service  adopted,  and  currently  uses,  a  method  that 
includes  the  work  done  by  Ferrel  (1886)  and  Bigelow  (1900).  The 
pressure  reduction  is  accomplished  using  the  hypsometric  equation  with 
the  temperature  replaced  by  the  mean  virtual  temperature  of  a 
fictitious  air  column  representing  conditions  between  the  station  and 
sea  level.  The  calculation  of  the  mean  virtual  temperature  includes 
computation  of  :  (1)  the  temperature  argument  defined  by  Bigelow 

(1900)  as  the  average  of  the  current  temperature  and  the  temperature 
recorded  twelve  hours  earlier,  (2)  a  standard  lapse  rate  correction, 

(3)  a  humidity  correction  first  introduced  by  Bigelow  (1900),  and  (4)  a 
plateau  correction  developed  by  Ferrel  (1886)  (Manual  of  Barometry, 
1963). 

Even  with  the  corrections  made  by  Ferrel  and  Bigelow,  the  National 
Weather  Service  pressure-reduction  method  still  includes  errors  that 
result  from  the  estimation  of  a  lapse  rate  through  the  plateau  or 
mountain. 

Bellamy  (1945)  attempted  to  mitigate  the  problem  of  estimating  a 
lapse  rate  in  developing  his  altimeter  correction  system  because  the 
altimeter  equation  assumes  a  standard  lapse  rate.  The  obvious  problem 
with  this  Is  that  the  atmosphere  seldom  follows  the  standard  lapse 
rate. 

To  sidestep  the  sea-level  reduction  difficulties,  Sangster  (1960) 
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devised  a  method  for  directly  computing  the  horizontal-pressure 
gradient  force  (HPGF)  at  the  earth's  surface  using  station  pressures 
and  surface  virtual  temperatures.  The  method  employs  Bellamy's 
altimeter  correction  system  in  order  to  rewrite  the  hydrostatic 
equation,  a  o-coordinate  system  developed  by  Phillips  (1957)  where  the 
vertical  coordinate  is  defined  by  pressure  with  the  surface  pressure 
used  as  the  lower  boundary  <o=l  at  the  surface),  a  smoothed  terrain  and 
finite  differences  to  produce  stream  and  potential  functions  of  the 
surface  geostrophic  wind  (Sangster,  1960). 

In  his  case  studies,  Sangster  found  the  geostrophic  wind  to  be 
less  when  computed  from  the  streamfunction  and  potential  function 
charts  than  when  computed  from  the  sea-level  charts  (Sangster,  1960). 

He  also  found  evidence  that  the  observed  winds  are  not  as  ageostrophic 
as  the  sea-level  isobars  would  indicate.  These  two  findings  result 
from  the  potential  function,  or  divergent  part  of  the  geostrophic  wind. 

In  isobarlc  coordinates,  this  divergent  part  of  the  wind  is  zero,  but 
since  the  sigma-coordinate  system  is  being  used,  there  is  a  divergence 
of  the  geostrophic  wind  and  thus  a  need  for  a  potential  function 
(Sangster,  1987).  Sangster's  case  studies  show  that  this  potential 
function  can  yield  a  component  in  almost  the  exact  opposite  direction 
as  that  given  by  the  streamfunction  in  some  locations,  thus 
illustrating  that  this  potential  function  is  not  negligible  (Sangster, 

1987). 

Another  method  for  deriving  a  pressure  fieid  and  subsequent 
geostrophic  wind  field  for  flat  terrain  was  developed  by  Pielke  and 
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Cram  (1987).  Using  a  terrain  following  coordinate  system  derived  by 
Pielke  and  Martin  (1981)  and  the  slope-flow  equations  developed  by 
Pielke  et  al  (1985),  Pielke  and  Cram  redefine  the  geostrophic  wind 
equations  where  the  pressure-gradient  term  is  the  sum  of  the  modified 
pressure  gradient  along  a  o-surface  and  terrain-gradient  terms.  This 
method  derives  a  flat-ground  surface-pressure  field  which  is  consistent 
in  concept  with  the  currently  applied  mean  sea  level  (MSL)  analyses, 
but  which  reduces  the  apparently  excessive  pressure  gradients  that 
result  from  conventional  MSL  procedures  (Pielke  and  Cram,  1987). 

Both  Sangster  and  Pielke  and  Cram  seem  to  have  circumvented  the 
difficulty  of  finding  a  temperature  for  a  fictitious  column  of  air. 

The  question  is  which  method  should  be  employed.  Davies-Jones  (1988) 
compares  both  methods  and  concludes  that  both  methods  are  equivalent. 
Davies-Jones  also  concludes  that  Sangster's  method  with  its  use  of  the 
altimeter  correction  system  and  introduction  of  a  horizontally  varying 
reference  atmosphere  is  more  complicated  than  the  method  devised  by 
Pielke  and  Cram  (Davies-Jones,  1988).  It  is  for  this  reason  that  this 
research  will  employ  the  Pielke  and  Cram  method  for  determining 
horizontal  pressure  gradient  over  flat  terrain. 

1.2  Research  Goals 

The  goals  of  this  research  are  to  use  the  data  collected  during 
the  Genesis  of  Atlantic  Lows  Experiment  (GALE)  to: 

(i)  Determine  corrections  to  the  pressure  data  recorded  at  the 
Portable  Automated  Mesonet  II  (PAM  II)  stations. 

(ii)  Determine  a  method  of  reducing  the  corrected  PAM  II  station 
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pressures  to  sea  level  In  order  to  enable  the  PAM  II  pressure  data  to 
be  used  in  mesoscale  analyses. 

(ill)  Study  methods  to  produce  a  meaningful  surface-pressure 
field  from  which  to  compute  surface  geostrophic  winds,  and  to  determine 
the  true  ageostrophic  components  of  the  winds  during  cold-air  damming 
situations. 

1.3  GALE  Project 

The  field  phase  of  the  GALE  project  was  conducted  from  15  January 
to  15  March  1986.  Some  of  the  objectives  of  GALE  (Dirks  et  al.,  1988) 
were  to: 

(i)  describe  the  airflow,  mass,  and  moisture  fields  in  East  Coast 
winter  storms  with  special  emphasis  on  mesoscale  processes. 

(li)  understand  the  physical  mechanisms  controlling  the  formation 
and  rapid  development  of  East  Coast  storms. 

(iii)  apply  conceptual  models  and  enhanced  mesoscale  data  to  the 
short  range  forecasting  of  significant  weather. 

To  accomplish  these  objectives,  analysis  of  mesoscale  pressure 
fields  Is  vital.  A  number  of  mesoscale  phenomena  and  mesoscale 
processes  are  associated  with  East  Coast  storms  and  their  roles  in 
cyclogenesis  needs  to  be  studied  (Dirks  et  al.,  1988).  Observations 
from  the  PAM  II  stations  and  near-shore  buoys  help  to  spatially  resolve 
borderline  meso-a  and  meso-p  scale  features.  The  6-minute  data 
recorded  at  the  PAM  II  sites  permit  the  detection  of  small  meso-p  and 
large  meso-r  scale  features  and  enable  their  observation  for  periods  of 
approximately  twelve  hours  (GALE,  1985). 
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2.  CASE  STUDY  DATA  AND  METHODOLOGY 

2.1  DATA 

The  data  for  this  research  are  a  combination  of  standard  hourly 
observations  at  National  Weather  Service  (NWS),  Federal  Aviation 
Administration  (FAA),  and  US  military  stations  and  special  observations 
at  selected  sites  operated  specifically  for  the  GALE  field  program. 

This  section  concentrates  on  the  special  observing  facilities  used  in 
this  research. 

The  objectives  of  GALE  require  meso-  and  even  microscale  analysis, 
so  the  data  collection  network  was  intensive.  There  are  two  data- 
gathering  areas  of  particular  significance  to  the  present  research: 

i)  Inner  GALE  Area 

In  this  area,  the  meso-|3  processes  of  frontogenesis  and 
cyclogenesis  can  be  examined.  The  inner  GALE  area  was  approximately 
500  km  wide,  centered  on  the  coast,  and  extended  1000  km  from  Georgia 
to  Virginia  (Figure  2.1).  Portable-Automated-Mesonet  (PAM  II),  Doppler 
radars,  ships,  buoys,  most  aircraft  flights  and  the  Cross-chain  Loran 
Atmospheric  Sounding  System  (CLASS)  rawinsonde  sites  were  located  in 
this  area. 

11)  Regional  GALE  Area 

The  regional  GALE  area  was  1000  km  wide  (from  the  Appalachian 
Mountains  to  600  km  offshore),  and  1500  km  long  (from  Florida  to  New 
Jersey;  see  Figure  2.1).  The  regional  area  will  allow  meso-a  processes 
of  frontogenesis  and  cyclogenesis  to  be  examined. 
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Figure  2.1  This  figure  shows  the  locations  of  the  Inner  GALE  Area 
(shaded  region)  and  the  Regional  GALE  Area. 
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2.1.1  Surface  Measurements 

The  surface  measurements  were  designed  to  provide  surface 
data  fields  of  standard  meteorological  parameters  within  the  Inner  GALE 
Area  with  mesoscale  resolution  and  to  provide  a  complement  to  the  data 
gathered  through  sounding  operations.  Standard  measurements  included 
air  temperature,  dew  point  temperature,  barometric  pressure,  wind  speed 
and  wind  direction.  In  addition,  land-based  stations  measured 
precipitation  and  sea-based  stations  measured  sea  temperature. 

The  fifty-station  PAM  II  network  (Figure  2.2)  provided 
automated  6-minute  averaged  meteorological  surface  observations  over 
the  eastern  half  of  North  Carolina,  all  of  South  Carolina,  and 
southeastern  Virginia.  A  line  of  four  stations  extended  northwestward 
to  support  a  sounding  cross-section  for  cold-wedge  studies. 

Deployment  of  the  eight  special  GALE  buoys,  including  six 
North  Carolina  State  University  buoys  and  two  NOAA-E  buoys  (Figure  2.2), 
supported  studies  of  the  development  of  the  coastal  front  and  augmented 
observations  in  the  data-sparse  region  in  the  Inner  GALE  Area. 

2.1.2  Ships 

The  oceanographic  research  vessels  R/V  Cape  Hatteras  (RVC) 
and  the  R/V  Endeavor  (RVE)  (the  R/V  Endeavor  participated  only  in 
January)  operated  off  the  coast  of  North  and  South  Carolina.  The  data 
recorded  Included  wind  speed  and  direction,  air  temperature,  sea- 
surface  temperature,  barometric  pressure,  and  humidity.  The  two  ships 
also  served  as  launch  platforms  for  soundings. 


Figure  2.2.  Map  showing  the  locations  of  the  PAM  II  stations  and 
marine  buoys,  ships  ,  and  platforms. 
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2.2  Barnes  Objective  Analysis  Scheme 

The  calculation  of  the  pressure  corrections  and  the  geostrophic 
wind  fields  for  this  study  is  aided  by  the  Barnes  (1964,  1973) 
objective  analysis  method.  The  Barnes  objective  analysis  scheme  is 
chosen  because  of  its  simplicity,  accuracy,  cost  effectiveness,  and  its 
several  advantages  over  other  similar  techniques.  From  Harms  (1985) 
the  advantages  of  the  Barnes  scheme  are: 

1)  The  scheme  is  computationally  simple  and  assigns  weights 
solely  as  a  function  of  distance  between  datum  and  grid  point. 

ii)  The  weight  parameter,  which  determines  the  degree  of 
smoothing,  can  be  chosen  prior  to  the  analysis  so  that  the  pattern 
resolvable  by  the  data  distribution  will  be  revealed. 

iii)  Since  the  weighting  function  approaches  zero  asymptotically 
with  distance,  the  Influence  of  data  may  be  extended  any  distance  to 
insure  that  a  sufficient  number  of  observations  Influence  each  grid 
point  without  changing  the  weight  function  and  ,  therefore,  the 
response  characteristics. 

iv)  Only  two  passes  through  the  data  are  needed  to  achieve  the 
desired  pattern. 

v)  "Noise"  is  adequately  filtered  from  the  analysis  so  further 
smoothing  by  application  of  additional  numerical  filters  is  not 
necessary. 

The  Barnes  objective  analysis  scheme  accepts  data  from  observation 
points  and  mathematically  interpolates  values  to  a  two-dimensional  grid 
array.  The  values  at  the  grid  points  can  be  used  to  interpolate  data 


to  station  locations  such  as  was  done  in  the  pressure  correction 
procedure  or  be  used  in  finite  difference  calculations  such  as  in 
calculating  the  geostrophic  wind  components.  An  additional 
modification  to  the  scheme  was  employed  in  determining  the  pressure 
corrections  enabling  the  program  to  interpolate  back  to  the  Individual 
PAM  II  locations.  For  example,  if  q  represents  any  meteorological 
variable,  the  Interpolated  value  is  Just  the  weighted  mean  (q)  of 
observations  surrounding  the  grid  point.  That  is, 

i=l 

q  =  -  (2.1) 

& 

i=i 


where  n  is  the  total  number  of  stations  Influencing  a  given  grid  point 
and  w  is  the  observation  weight. 

The  observation  weights  are  assigned  to  the  data  as  a  known 
function  of  distance  between  the  data  and  grid  points  according  to 

w  =  exp(-dVk)  (2.2) 

where  k  la  the  weight  parameter  that  controls  the  rate  at  which  the 
weight  value  decreases  outward  from  the  point  of  interpolation.  Hence, 
k  determines  the  degree  of  smoothing  of  the  data  field.  If  k  is  small, 
there  is  little  smoothing;  If  k  is  large,  there  is  greater  smoothing. 

The  selection  of  k  is  important  in  determining  the  detail 
remaining  in  the  Interpolated  field.  The  choice  of  this  parameter 
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value  must  strike  a  balance  between  retaining  as  much  detail  as  the 
observation  network-density  allows  and  filtering  out  sources  of  random 
and  sub-scale  error.  Structural  detail  Is  limited  by  the  minimum 
resolvable  wavelength  which  is  usually  twice  the  station  spacing.  The 
GALE  PAM  II  stations  have  a  mean  station  spacing  of  about  68  km,  while 
the  upper  air  stations  of  the  Inner  GALE  sounding  (including  National 
Weather  Service  and  CLASS)  network  have  an  average  station  spacing  of 
about  180  km.  These  distances  are  about  one-half  those  of  the  normal 
operational  reporting  networks.  Thus,  theoretically,  the  Inner  GALE 
networks  can  resolve  features  of  twice  the  respective  mean  station 
separations,  or  about  140  km  at  the  surface  and  360  km  aloft. 

The  weight  parameter,  k.  is  selected  to  determine  the  percentage 
of  the  amplitudes  of  the  minimally-resolved  waves  that  will  be 
conserved  (Harms,  1985).  The  specific  k  values  used  in  this 
application  suppress  the  amplitudes  of  the  two-station -spacing  waves  to 
only  16  percent  of  their  original  amplitude  in  the  belief  that  error 
does  make  a  significant  contribution  to  the  amplitudes  of  the  higher 
frequency  features.  Sources  of  "error"  Include  biases  introduced  by 
local  topography  and  obstacles  at  observation  sites,  features  on  scales 
smaller  than  the  station  spacing  (for  example,  gravity  waves),  and 
turbulent  fluctuations  with  periods  of  several  minutes  (especially  near 
the  ground).  Specific  applications  of  the  Barnes  analysis  scheme  will 
be  described  in  later  sections. 
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3.  PRESSURE  CORRECTION  AND  REDUCTION  TO  SEA  LEVEL 

The  data  recorded  at  the  PAM  II  sites  will  enable  researchers  to 
gain  a  better  understanding  of  the  mesoscale  processes  involved  in  East 
Coast  winter  storms.  A  problem  with  the  PAM  II  data  is  that  the 
recorded  pressure  has  not  been  corrected  or  reduced  to  a  horizontal 
reference  level  and,  as  a  result,  cannot  be  used  directly  to  perform 
mesoscale  analyses.  This  section  first  presents  one  method  of  finding 
pressure  corrections  for  each  PAM  II  station  and  then  describes  a 
method  for  reducing  these  corrected  pressures  to  sea  level. 

3.1  Pressure  Correction 

Correction  of  PAM  II  station  pressures  involves  a  series  of  steps. 
First,  buddy  checks  were  performed  between  altimeter  settings  at 
selected  PAM  II  sites  and  nearby  NWS  stations.  Second,  time  series 
plots  of  altimeter  differences  for  each  selected  station  pair  were 
evaluated  in  order  to  identify  trends.  Third,  and  finally,  an 
objective  analysis  based  only  on  the  NWS  altimeter  settings  was 
employed  to  determine  the  corrections  for  all  PAM  II  stations. 

Initially  five  comparisons  between  the  calculated  altimeter 
settings  of  PAM  II  and  NWS  stations  were  made  for  24  January,  24-25 
February,  and  8  March  using  PAM  II  sites  3,  9,  26,  33,  34  and  NWS 
stations  AGS,  CAE,  GSO,  and  RDU.  No  long-term  time  trends  in  the 
differences  were  found  for  these  three  periods,  so  the  study  was 
expanded  to  include  the  entire  GALE  period. 

Buddy  checks  were  performed  for  the  entire  GALE  period  between 
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each  of  ten  PAM  II  stations  and  its  nearest  NWS  neighbor  using  a 
calculated  altimeter  setting  for  each  station.  The  altimeter  setting 
was  chosen  in  order  to  reduce  the  pressure  offsets  due  to  elevation 
differences  between  the  station  pairs.  The  equation  used  to  calculate 
the  altimeter  setting  is  taken  from  the  Smithsonian  Meteorological 
Tables  (List,  1966)  and  is  shown  below: 

Alstglb=(pik-0.3)1 1  +<ap01/T0)(H/(pirO.3)*)l1/\  (3.1) 

where:  Alstg^  =  altimeter  setting  in  mb, 
plb  =  station  pressure  in  mb, 
p„  =  standard  sea  level  pressure  =  1013.25  mb, 
a  =  lapse  rate  in  standard  atmosphere  =  0.0065’C/m, 

T0  =  standard  sea  level  temperature  =  288  K, 

H  =  station  elevation  in  meters, 

n  =  aR  =  0.190284,  where  R  is  the  dry-air  gas  constant. 

Table  3.1  shows  the  results  of  the  station  comparisons  over  the 
entire  GALE  period. 

In  order  to  determine  how  much  accuracy  could  be  expected  for  any 
given  PAM  II  barometer,  the  NCAR  field  calibration  checks  for  the 
pressure  sensors  were  examined.  During  field  calibrations,  teams 
transported  a  barometer  calibrated  at  the  nearest  National  Weather 
Station  to  a  nearby  PAM  II  site  and  compared  the  two  instruments.  Each 
PAM  II  barometer  was  checked  at  least  two  and  at  most  five  times  during 
GALE.  Most  of  the  sites  had  at  least  three  field  calibration  checks 
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Table  3.1.  Station  Comparisons  with  the  distance  between  stations  in 
km,  and  average  altimeter  difference,  standard  deviation,  and  range  in 
millibars. 


STATION  COMPARISONS 


DISTANCE 

DIFFERENCE 

STD  DEV 

RANGE 

NWS 

PAM 

(KM) 

(MB) 

(MB) 

(MB) 

SAV 

01 

25 

0.281 

0.378 

-1.02  TO 

1.41 

AGS 

03 

57 

-0.837 

0.686 

-2.95  TO 

1.03 

CHS 

06 

35 

0.544 

0.599 

-2.54  TO 

2.86 

CAE 

09 

20 

-1.340 

0.458 

-2.51  TO 

0.10 

CLT 

15 

54 

-1.045 

0.726 

-0.72  TO 

2.48 

ILM 

21 

50 

2.309 

0.705 

0.61  TO 

4.33 

GSO 

26 

25 

-0.649 

1.003 

-4.08  TO 

2.25 

RDU 

33 

29 

1.763 

0.448 

0.48  TO 

3.35 

RDU 

34 

33 

1.188 

1.194 

-0.51  TO 

2.81 

GSO 

34 

75 

-0.098 

1.355 

-3.67  TO 

2.81 

HAT 

45 

13 

-0.907 

0.475 

-4.39  TO 

0.60 

performed  during  the  GALE  period. 

The  field  calibration  data  were  not  directly  used  to  correct  the 
PAM  II  station  pressures  because,  as  Wade  (1985)  states,  the  direct 
comparison  with  reference  sensors  provides  a  check  of  the  absolute 
accuracy  of  the  mesonet  sensors,  but  the  detection  of  errors  in  the 
data  involves  the  comparison  with  similar  sensors  in  the  network.  In 
other  ^ords,  in  order  to  fit  the  mesonet  into  the  existing  NWS  network, 
the  mesonet  sensors  must  be  compared  with  the  surrounding  NWS  stations 
to  check  for  errors. 

Since  it  is  difficult  to  obtain  a  reliable  statistic  for  any  given 
site  based  on  the  few  comparisons  performed  per  station,  all  of  the 
comparisons  for  all  sites  were  combined  as  follows:  An  average  field 
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calibration  difference  was  calculated  for  each  site.  Then,  using  this 
average,  an  absolute  deviation  from  the  mean  difference  was  calculated 
for  each  calibration  time  for  each  station.  These  deviations  were 
averaged  over  the  entire  PAM  II  network  to  obtain  an  average  deviation 
for  the  entire  network  of  0.375  mb,  which  represents  a  measure  of  the 
absolute  resolution  of  the  PAM  II  sensors.  Comparing  this  with  the 
average  differences  in  Table  3.1  indicates  that  most  of  the  stations 
can  be  improved. 

The  next  step  was  to  plot  the  altimeter  differences  over  time  to 
evaluate  the  trends.  Over  the  sixty-day  period,  the  differences  show 
variations  of  +/-  0.5  to  1.0  mb  over  the  synoptic  2  to  7  day  period, 
and  occasional  sharp  +/-  1  to  2  mb  variations  over  2-  to  3-  hour 
periods. 

Comparison  with  daily  sea-level  pressure  analyses  shows  that, 
generally,  the  differences  remain  fairly  well-behaved  during  high 
pressure  situations  while  extended  deviations  from  the  mean  occur 
during  periods  when  low-pressure  systems  are  in  the  area.  The 
magnitude  of  the  differences  depends  on  the  orientation  of  the  PAM/NWS 
pair  with  the  pressure  gradient.  For  example,  if  the  pressure  gradient 
was  oriented  north-south,  two  stations  oriented  east-west  would  have 
smaller  differences  than  two  stations  oriented  north-south. 

Figure  3.1  shows  an  example  of  the  synoptic  variations.  In  this 
figure,  the  differences  fall  below  the  mean  (0.281  mb)  around  0000  UTC 
on  19  January  and  remain  below  the  mean  until  around  1200  UTC  on  21 
January.  During  this  time,  a  low  pressure  center  formed  north  of  the 
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Figure  3.1  Plot  of  the  altimeter  difference  between  Savannah,  GA  (SAV) 
and  PAM  II  Station  #01  with  time  the  period  0600  UTC  on  18  January 
1986  to  1800  UTC  on  21  Janupry  1986. 

Inner  GALE  region  and  a  persistent  surface  trough  was  located  west  of 
SAV. 

Attention  is  now  turned  to  the  short  period  variations  to 
determine  if  they  could  be  eliminated  from  the  data  set.  Figure  3.2  Is 
a  good  example  of  a  short  period  deviation.  The  differences  sharply 
decrease  from  1600  UTC  6  February  to  a  low  of  -1.60  mb  at  1800  UTC  on  6 
February  and  return  to  the  mean  two  hours  later.  Figure  3.3,  showing 
the  altimeter  settings  for  both  SAV  and  PAM  II  station  #01,  helps  to 
explain  this  spike.  The  sharp  decrease  In  the  difference  is  a  result 
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Figure  3.2  Same  as  Figure  3.1  except  for  the  period  2200  UTC  on  5 
February  1986  to  1400  UTC  on  7  February  1986. 

of  a  temporary  rise  in  the  altimeter  setting  at  SAV  while  the  altimeter 
setting  at  the  PAM  II  station  continued  to  drop.  A  frontal  system 
approaching  from  the  west  would  explain  the  continued  drop  in  altimeter 
setting,  but  a  check  of  the  observations  at  SAV  reveals  that  a 
thunderstorm  was  observed  at  1610  UTC,  lasting  for  about  forty  minutes. 
The  resulting  roesohigh  lasted  for  two  observations  (1700  UTC  and  1800 
UTC)  before  moving  off  and  allowing  the  pressure  to  fall  rapidly  as  the 
frontal  system  approached.  Therefore,  the  spike  depicted  in  Figure  3.2 
is  likely  the  result  of  a  real  small-scale  feature  and  thus  it  can  not 
be  eliminated  from  the  data  set.  The  other  spikes  in  the  altimeter 
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Figure  3.3.  Plot  of  the  Altimeter  Settings  for  SAV  and  PAM  II  station 
#01  from  2200  UTC  on  5  February  1986  to  1400  UTC  on  7  February  1986. 


difference  plot  are,  as  in  Figure  3.2,  not  usually  composed  of  isolated 
points  so  it  was  decided  to  retain  all  of  the  field  data. 

Since  not  all  PAM  II  sites  are  located  near  NWS  stations,  the 
final  step  in  the  correction  process  involves  performing  an  objective 
analysis  on  a  field  of  NWS,  FAA,  military,  and  marine  altimeter 
settings  to  obtain  corrections  for  each  station. 

A  map  showing  the  PAM  II  sites  and  the  NWS,  FAA,  military,  and 
marine  stations  used  in  the  Barnes  objective  analysis  is  shown  in 
Figure  3.4.  The  twelve  marine  observations  help  to  define  the 
pressure  field  over  the  nearby  waters.  Data  from  nearshore  buoys  were 
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eliminated  because,  during  the  data  reduction  and  calibration  phase  of 
GALE,  the  NCSU  buoys  were  found  to  not  be  as  reliable  as  expected  due 
to  unpredictable  drift  in  calibration.  This  elimination  left  only  the 
twelve  stations  shown  in  Figure  3.4.  Since  the  height  of  these  marine 
observations  is  taken  as  zero,  Equation  (3.1)  reduces  to: 

Alstgtk  =  (p>k  -  0.3).  (3.2) 

The  cases  chosen  for  analysis  have  a  variety  of  pressure  gradients 
in  order  to  obtain  a  better  representation  of  the  entire  GALE  period. 

Five  cases  (24  January,  14  and  23  February,  and  1  and  8  March)  are 
characterized  by  small  gradients,  three  cases  (2.  13,  and  25  February) 
have  a  moderate  gradient,  two  (20  January  and  6  March)  have  a  strong 
gradient,  and  three  cases  (11  and  27  February  and  13  March)  have  low 
pressure  systems  in  the  GALE  region.  Thirteen  cases  were  selected  so 
that  there  would  be  at  least  nine  offset  values  for  each  PAM  II  site 
which  would  then  be  averaged  to  obtain  the  correction  for  each  station. 

The  1200  UTC  hour  was  chosen  to  be  used  in  each  case  because  it  is 
a  synoptic  hour  when  a  sufficient  amount  of  NWS  and  marine  data  are 
available  to  perform  the  Barnes  analysis.  Furthermore,  in  order  to 
check  the  consistency  of  the  Barnes  analysis  program,  the  program  was 
run  on  data  for  five  consecutive  hours  starting  at  1000  UTC  on  14 
February. 

A  check  was  performed  before  applying  the  Barnes  analysis  scheme 
In  order  to  Identify  any  bogus  data.  The  altimeter  settings  of  each 
station  were  plotted  for  each  case.  Upon  examination  of  the  resulting 


Figure  3.4.  Map  showing  the  surface  stations  employed  in  the  Barnes 
Objective  Analysis  Scheme.  The  PAM  II  stations  are  included  as  a 
reference. 
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altimeter  field,  the  altimeter  settings  at  the  military  stations  NBC, 

NCA,  and  NKT  were  found  to  be  consistently  lower  than  the  surrounding 
stations.  As  a  result,  these  three  stations  were  removed  from  the  data 
set  used  by  the  Barnes  analysis  scheme. 

The  Barnes  analysis  scheme  was  then  applied  to  the  remaining 
surface  NWS,  FAA,  military,  and  marine  stations.  The  mean  station 
spacing  for  this  network  is  about  130  km.  Thus,  it  follows  that  the 
minimum  resolvable  waves  are  those  of  260  km  wavelength.  A 
corresponding  weight  parameter,  k,  of  10*  km1  is  used  to  suppress  these 
minimum  resolvable  waves  to  about  16  percent  of  their  initial  amplitude 
at  the  end  of  the  second  pass  of  the  scheme.  In  order  to  ensure 
convergence  of  the  scheme  after  two  passes,  a  convergence  parameter  (D 
of  0.3  is  used.  Altimeter  settings  are  interpolated  to  a  two- 
dimensional  array  of  grid  points  with  a  mean  station  spacing  of  40  km 
at  60*N  or  about  34  km  at  the  PAM  II  network  latitude. 

After  the  second  pass  of  the  Barnes  scheme,  an  altimeter  setting 
resulted  at  each  grid  point  as  shown  in  Figure  3.5.  In  order  to 
determine  the  altimeter  setting  at  each  PAM  II  station,  another 
interpolation  was  done  using  the  data  from  the  four  surrounding 
grid  points.  This  interpolation  was  accomplished  by  bilinear 
interpolation  from  the  surrounding  grid  points. 

Once  an  altimeter  setting  was  obtained  for  each  PAM  II  site,  the 
program  compared  this  altimeter  setting  with  an  altimeter  setting 
calculated  from  the  PAM  II  height  and  raw  station  pressure  data  using 
equation  (3.1).  An  offset  was  calculated  between  the  interpolated  and 
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analysis  scheme. 

observed  values  by  using  the  following  equation: 

DIFF  =  ALSTGob  -  ALSTG^.  (3.3) 

The  results  for  the  thirteen  cases  are  shown  in  Table  3.2  and  the 
results  of  the  five  consecutive  hours  are  shown  in  Table  3.3. 

From  Table  3.3,  the  corrections  obtained  from  the  Barnes  scheme 
for  nearly  all  stations  during  the  five  consecutive  hours  range  within 
+/-  0.375  mb,  which  was  previously  determined  from  the  field 
calibration  data  to  be  the  resolution  of  the  PAM  II  instruments. 
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Table  3.2.  This  table  shows  the  results  of  the  Barnes  objective 
analysis.  Units  for  the  average  differences  and  standard  deviation  of 
the  differences  are  in  aillibars. 


Station 

DIFF 

c 

Station 

DIFF 

0 

01 

0.274 

0.216 

26 

0J46 

0.769 

02 

0.566 

0.336 

27 

0.712 

0.376 

03 

-1.195 

0.293 

28 

-1.509 

0.788 

04 

-0.055 

0.349 

29 

-1.363 

0.533 

05 

1.075 

0.413 

30 

-0.069 

0.443 

06 

0.569 

0.681 

31 

-0.040 

0.468 

07 

1.361 

0.473 

32 

2.505 

0.600 

08 

0.712 

0.335 

33 

0.146 

0.267 

09 

-0.880 

0.422 

34 

-0.220 

0.249 

10 

-1.833 

0.727 

35 

2.461 

0.528 

11 

-0.120 

0.470 

36 

0.267 

0.386 

12 

0.261 

0.566 

37 

0.072 

0.335 

13 

0.718 

0.353 

38 

-0.072 

0.436 

14 

-0.091 

0.259 

39 

-0.513 

0.325 

15 

-0.318 

0.344 

40 

-0.145 

0.285 

16 

-0.577 

0.276 

41 

1.899 

0.356 

17 

0.335 

0.349 

42 

0.162 

0.222 

18 

0.171 

0.319 

43 

-0.618 

0.508 

19 

-1.025 

0.197 

44 

0.156 

0.597 

20 

1.171 

0.461 

45 

-0.804 

0.449 

21 

2.088 

0.437 

46 

0.070 

0.432 

22 

1.568 

0.216 

47 

0.035 

0.234 

23 

-0.087 

0.349 

48 

-0.040 

0.295 

24 

1.538 

0.965 

49 

-0.739 

0.330 

25 

0.952 

0.360 

50 

-0.268 

0.486 

Therefore,  it  would  appear  that  the  Barnes  objective  analysis  scheme  Is 
producing  valid  results. 

Table  3.4  shows  a  comparison  between  the  ten  station  buddy  check 
results  and  the  average  differences  obtained  with  the  Barnes  analysis 
for  the  same  ten  stations.  All  but  two  of  the  stations  have  the  same 
sign.  An  explanation  for  the  change  In  sign  for  PAM  #26  is  that  the 
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Title  3.3  Tkii  title  tiois  tie  correction  oktiiied  froi  the  lines  uilysii  iclne  for  tie  file  eoueeitiie  loirs 
kefiuiif  lit!  1110  fffC  oi  14  Felrnry  1)34.  Blits  for  tie  correction  ire  lillilirs. 


ftitioi 

1000 

1100 

1200 

1300 

1400 

Stitici 

1000 

1100 

1200 

1300 

1400 

(1 

-0.11 

-0.21 

0.14 

0.01 

0.14 

24 

-0.2) 

-0.14 

-0.35 

-0.11 

0.0) 

02 

0.3) 

0.53 

0.)1 

0.21 

0.44 

27 

1.04 

0.)5 

1.32 

1.75 

1.54 

03 

-1.23 

-1.01 

-0.)5 

-1.2) 

-0.)4 

21 

-1.)) 

-2.14 

-1.43 

-1.72 

-l.)0 

14 

-0.3) 

0.15 

-0.17 

-0.23 

-0.20 

2) 

-2.13 

-2.04 

-1.54 

-1.15 

-1.5) 

15 

1.20 

0.)0 

0.73 

0.40 

0.41 

30 

-0.37 

-0.15 

-0.54 

0.31 

0.24 

04 

0.45 

0.05 

0.43 

0.25 

0.53 

31 

0.44 

0.03 

0.44 

-0.43 

0 

07 

0.)2 

1.2) 

1.15 

0.4) 

1.25 

32 

2.11 

1.24 

2.01 

l.)5 

1.13 

01 

0.4) 

0.43 

1.01 

0.42 

0.45 

33 

0.20 

-0.14 

0.32 

0.10 

0.14 

05 

-1.22 

-1.07 

-1.24 

-1.43 

-1.34 

34 

-0.7) 

-0.42 

-0.51 

-0.34 

-0.30 

10 

-1.42 

-1.34 

-1.52 

-1.72 

-1.2) 

35 

1.7) 

l.)7 

1.54 

2.41 

2.2) 

11 

-0.24 

-0.)5 

-0.24 

-0.47 

-0.47 

34 

-0.04 

0.03 

-0.10 

0.34 

0.27 

12 

-0.24 

0.45 

0.52 

0.3) 

0.24 

17 

-0.20 

-0.23 

-0.7) 

0.03 

-0.21 

13 

0.54 

1.04 

1.17 

0.43 

0.40 

34 

-0.31 

-0,)2 

0.21 

0.11 

-0.04 

14 

0.33 

0.04 

0.04 

-0.01 

0.1) 

3) 

1.02 

0.)6 

-0.2) 

-1.04 

-1.34 

15 

-0.37 

-0.44 

-0.2) 

-0.33 

-0.07 

40 

-0.71 

-0.74 

-0.22 

-0.45 

-0.25 

14 

-1.47 

-1.0) 

-0.25 

-0.44 

-0.7) 

41 

1.43 

1.44 

1.44 

1.44 

1.73 

IT 

-0.14 

0.13 

0.10 

0.02 

-0.12 

42 

-0.44 

-0.43 

-0.17 

-0.07 

-0.31 

11 

-0.45 

0.35 

0.0) 

0.20 

-0.42 

41 

-0.55 

-0.75 

-0.47 

-1.00 

-1.10 

11 

-1.05 

-0.44 

-1.23 

-1.4) 

-1.30 

44 

-0.44 

0.20 

0.33 

-0.0) 

0.01 

20 

0.)4 

0.)) 

1.15 

1.43 

1.17 

45 

-1.34 

-l.)4 

-1.44 

-1.74 

-0.59 

21 

2.22 

2.02 

2.30 

2.15 

1.J7 

44 

-0.34 

-0.55 

-0.33 

-0.51 

0.37 

22 

0.14 

0.02 

1.23 

0.74 

1.0) 

47 

0.21 

-0.47 

-0.10 

0.33 

0.71 

23 

-0.27 

-0.13 

-0.4) 

-0.24 

-0.27 

41 

-0.05 

-0.17 

-0.11 

-0.07 

0.31 

24 

liiit 

mtt 

1.5) 

1.44 

2.00 

4) 

-0.)5 

-0.44 

-0,73 

-o.n 

-0.41 

25 

0.45 

0.3) 

0.2) 

0.24 

1.04 

50 

-0.34 

-0.14 

-0.30 

-0.21 

-0.44 

station  is  closer  to  Winston-Salem  (INT)  than  Greensboro  (GSO)  so  the 
Barnes  analysis  would  weigh  more  heavily  the  Winston-Salem  altimeter 
setting  than  that  of  Greensboro.  Also,  one  would  expect  the  average 
differences  of  the  Barnes  analysis  and  the  buddy  check  to  differ 
because  the  former  method  incorporates  data  from  surrounding  stations 
into  the  interpolated  altimeter  setting  instead  of  data  from  only  one 
station.  Thus,  at  least  the  synoptic  and  meso-a  scale  pressure  field 
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Table  3.4  Comparison  between  the  differences  obtained  from  the  Buddy 
Check  and  Barnes  Analysis  methods.  Units  of  the  "corrections"  are  in 
millibars. 


Station 

Buddy  Check 

Barnes  Analysis 

01 

0.211 

0.274 

03 

-0.682 

-1.195 

06 

0.693 

0.569 

09 

-1.475 

-0.880 

15 

1.501 

-0.318 

21 

2.522 

2.088 

26 

-0.293 

0.148 

33 

1.730 

0.146 

34 

0.741 

-0.220 

45 

0.777 

-0.804 

EX 

5.725 

-0.192 

Mean 

0.573 

-0.019 

a 

1.128 

0.883 

are  accounted  for. 

One  might  expect  the  Barnes  method  which  accounts  for  the 
synoptic-scale  pressure  patterns  to  represent  an  improvement  over  the 
buddy  check  method.  In  order  to  determine  whether  there  is  a 
statistically  significant  difference  between  the  two  methods,  a  test  of 
the  variances  of  each  method  was  conducted.  A  Chi-squared  test  was 
performed  using  the  following  information: 

oK  =  1.128  mb  (standard  deviation  for  buddy  check  method) 
aik  ~  0.883  mb  (standard  deviation  for  Barnes  method) 

Level  of  significance,  a  =  0.05 

The  null  hypothesis  and  the  alternate  hypothesis  were  taken  as 
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H0:  °*BA  =  and 

Hj:  0>u  <  o*K 

From  Table  VI  in  Walpole  and  Meyers  (1978),  and  using  v  =  9 
degrees  of  freedom,  the  critical  region  was  taken  as  X*  <  3.325,  where 
X*  =  (n-l)o*gA/o*|Kj  In  this  test  if  a  calculated  X*  is  less  than 
3.326,  the  null  hypothesis  can  be  rejected  and  the  alternate  hypothesis 
accepted.  The  value  of  X*  was  found  to  be  5.515.  Since  this  is 
greater  than  the  critical  value  of  3.325,  the  null  hypothesis  can  be 
accepted  and  it  can  be  concluded  that  there  is  no  significant 
difference,  statistically,  between  the  two  methods.  However,  the 
Barnes  analysis  method  will  be  used  for  the  pressure  corrections 
because  not  all  of  the  PAM  II  stations  have  a  nearby  NWS  station  with 
which  to  cond  ,<  t  a  buddy  check. 

As  can  be  seen  from  Table  3.2,  some  of  the  differences  have 
absolute  values  greater  than  1  mb.  In  order  to  determine  if  these 
large  differences  are  the  result  of  PAM  II  elevation  error  or  Just 
instrument  offsets,  a  check  of  the  NCAR  field  calibration  data  was 
performed.  The  field  calibration  differences  were  averaged  for  each 
station  and  these  averages  were  compared  to  the  stations  whose  average 
Barnes  differences  exceed  +/-  1  mb.  This  comparison  is  presented  in 
Table  3.6.  In  most  cases  where  the  magnitude  of  the  Barnes  scheme 
correction  exceeds  1  mb,  the  average  field -calibration  difference  also 
exceeds  1  mb.  However,  some  stations  experience  a  sign  change  between 
the  field  calibrations  and  the  Barnes  analysis.  As  a  result,  this 
check  would  appear  to  be  inconclusive  on  whether  the  large  corrections 
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Table  3.5  Comparison  between  the  field  calibration  data  and  the  Barnes 
Analysis  data  for  selected  PAM  II  sites.  Units  are  millibars. 


Station 

Calibration 

Barnes 

03 

1.255 

-1.195 

05 

1.470 

1.075 

07 

2.140 

1.361 

10 

1.640 

-1.833 

19 

0.717 

-1.025 

20 

1.667 

1.171 

21 

2.473 

2.088 

22 

2.313 

1.568 

24 

2.506 

1.538 

28 

0.498 

-1.509 

29 

0.573 

-1.363 

32 

2.225 

2.505 

35 

2.793 

2.461 

41 

2.190 

1.899 

are  due  to  elevation  errors. 

A  second  check  into  this  matter  of  large  corrections  was  conducted 
by  examining  the  local  relief  of  those  stations  listed  in  Table  3.5. 

The  effect  of  local  relief  hinges  on  the  uncertainty  of  the  exact 
station  location  and  exact  terrain  elevation.  Once  the  station 
locations  were  found  on  a  topographic  map,  a  half-kilometer  diameter 
circle  was  drawn,  centered  on  the  station.  The  resolution  of  the 
topographic  maps  used  would  allow  the  station  location  to  be  pinpointed 
to  within  a  half-kilometer.  The  maximum  and  minimum  elevations  within 
this  circle  were  recorded  and  the  elevation  differences  were  found  and 
are  displayed  In  Table  3.6. 

From  Table  3.6,  only  stations  5,  19,  20,  28,  and  29  have  elevation 
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Table  3.6.  Shows  the  elevation  differences  (in  meters)  within  a  half¬ 
kilometer  circle  centered  at  the  station. 

Station  Difference 


03 

6 

05 

13 

19 

10 

20 

10 

21 

3 

22 

3 

24 

6 

28 

20 

29 

25 

32 

3 

35 

3 

41 

3 

differences  large  enough  to  account  for  corrections  exceeding  1  mb. 
These  five  stations  are  located  in  the  western  Piedmont  and  the 
Appalachian  regions  of  North  Carolina,  South  Carolina,  and  Virginia. 
Those  stations  located  in  the  Coastal  Plain  only  had  elevation 
differences  of  3-6  m,  which  would  not  account  for  their  large 
corrections.  From  this  it  would  appear  that  only  the  corrections  for 
the  stations  in  the  higher  elevations  may  be  the  result  of  errors  in 
station  elevation.  In  light  of  the  mixed  results,  it  seems  easier  to 
apply  the  large  corrections  to  the  station  pressure  rather  than  to 
correct  the  elevations  of  an  uncertain  number  of  PAM  II  stations. 

After  applying  the  correction  to  the  observed  altimeter  setting,  a 
new  station  pressure  was  calculated  using  the  following  equation: 


j 


P  =  |  (ALSTG)”  -  (P0”aH/T0)  )1/n  +  0.3.  (3.4) 


The  observed  station  pressure  was  subtracted  from  this  "corrected" 
station  pressure  to  find  the  pressure  correction  value.  Table  3.7 
gives  the  station  pressure  correction  for  each  PAM  II  station. 

3.2  Reduction  of  Station  Pressures  to  Sea  Level 

The  procedure  for  reducing  the  corrected  station  pressures  to  sea 
level  was  taken  directly  from  the  Manual  of  Barometry  (1963).  The 
procedure,  as  outlined  in  the  Manual  of  Barometry,  is  based  upon  the 
hypsometric  equation: 

P0  =  P  •  10(“™,  (3.5) 

where  P0  =  sea  level  pressure  (mb) 

P  =  station  pressure  (mb) 

K  =  hypsometric  constant  (0.0266895*R/m) 

H  =  station  elevation  (m) 

T>y  =  mean  virtual  temperature  CR) 

Use  of  the  hypsometric  equation  assumes  a  vertical  air  column  in 
which  the  effects  of  vertical  acceleration  are  negligible.  This 
assumption  describes  a  hydrostatic  condition  in  which  the  weight  of 
each  layer  In  the  air  column  is  in  balance  with  the  pressure  difference 
between  the  lower  and  upper  level  surfaces  confining  the  layer. 

The  Manual  of  Barometry  provides  tables  of  values  of  the  ratio 
P0/P=  r  which  are  based  on  the  elevation  and  mean  virtual  temperature 
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Table  3.7  Pressure  corrections  in  millibars  to  be  applied  to  the 
station  pressure  recorded  at  the  PAM  II  stations. 


Station 

Correction 

Station 

Correction 

01 

-0.27 

26 

-0.15 

02 

-0.57 

27 

-0.70 

03 

+1.18 

28 

+1.40 

04 

+0.05 

29 

+1.26 

05 

-1.07 

30 

+0.07 

06 

-0.57 

31 

+0.04 

07 

-1.36 

32 

-2.50 

08 

-0.71 

33 

-0.15 

09 

+0.86 

34 

+0.21 

10 

+1.80 

35 

-2.46 

11 

+0.12 

36 

-0.27 

12 

-0.26 

37 

-0.08 

13 

-0.71 

38 

+0.06 

14 

+0.08 

39 

+0.50 

15 

+0.31 

40 

+0.15 

16 

+0.58 

41 

-1.90 

17 

-0.34 

42 

-0.16 

18 

-0.17 

43 

+0.62 

19 

+1.01 

44 

-0.15 

20 

-1.16 

45 

+0.81 

21 

-2.09 

46 

-0.07 

22 

-1.57 

47 

-0.04 

23 

+0.09 

48 

+0.04 

24 

-1.57 

49 

+0.74 

25 

-0.94 

50 

+0.26 

of  the  station.  In  this  research,  r  was  calculated  using  the  equation: 

r  =  I0(n/™  (3.6) 

Virtual  temperature  is  used  in  this  procedure  in  order  to  account 
for  the  effect  of  water  vapor  on  the  density  of  moist  air.  According 
to  the  Manual  of  Barometry,  the  mean  virtual  temperature  can  be  found 
from  the  following  equation: 


(3.7) 
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T|?  =  T0  +  li  +  aH/2  +  e,Ch  +  F. 


where  Tfl  is  the  conversion  of  temperature  into  absolute  units,  either 
273  K  or  459. 7*R  and  t,  is  the  station  temperature  argument,  which  is 
simply  the  average  of  the  observed  temperature  and  the  temperature 
recorded  twelve  hours  earlier: 


t,  =  (t0  +  t.u)/2.  (3.8) 

The  standard  lapse  rate  correction  is  aH/2,  where  a  is  the 
standard  lapse  rate;  either  0.0065*C/m  or  0.117*F/m.  Using  this 
correction  means  that  the  mean  temperature  of  the  air  column  is 
considered  greater  than  the  temperature  argument  at  the  station  (t,)  by 
aH/2.  This  is  the  same  as  saying  the  mean  temperature  is  the  linear 
average  of  t,  and  the  temperature  at  sea  level,  based  on  the  assumption 
that  the  lapse  rate  has  a  standard  value. 

The  humidity  correction  is  e,^  with  e,  being  the  saturated 
vapor  pressure  (in  mb)  at  the  station  at  the  time  of  the  observation 
and  Ck  being  a  known  function  of  station  elevation  in  *F/mb.  This 
humidity  correction  is  necessary  because  the  density  of  moist  air  is 
less  than  the  density  of  dry  air  provided  that  the  pressure  and 
temperature  are  the  same.  The  quantity  e,Cj,  is  intended  to  represent  a 
correction  considered  as  a  difference  between  the  mear  virtual 
temperature  of  the  fictitious  air  column  and  the  mean  dry-bulb 
temperature  of  the  fictitious  air  column.  The  quantity  ef  refers  to 
the  surface,  but  the  quantity  e,^  refers  to  the  entire  air  column. 

The  correction  for  the  plateau  effect  and  local  lapse  rate  anomaly 
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is  F.  The  plateau  effect  correction  is  introduced  so  the  amplitude  of 
the  annual  variation  of  pressure  reduced  to  sea  level  is  approximately 
the  same  at  all  stations  regardless  of  elevation.  The  local  lapse  rate 
anomaly  correction  takes  into  account  any  deviation  of  the  assumed 
local  lapse  rate  from  the  standard  lapse  rate.  For  stations  with 
elevations  of  305  m  or  lower,  F  is  determined  on  the  basis  of  the 
annual  mean  station  temperature.  For  stations  exceeding  305  m  in 
elevation,  F  is  determined  by  geographic  location,  elevation,  and 
climatological  conditions. 

In  calculating  the  sea  level  pressure  for  each  PAM  II  station, 
estimates  of  the  annual  mean  temperature  and  saturated  vapor  pressure 
had  to  be  made.  An  annual  mean  temperature  for  each  station  was 
obtained  by  using  the  annual  mean  temperatures  of  the  nearby  NWS 
stations. 

The  saturated  vapor  pressure  was  obtained  from  Table  7.2.1  in  the 
Manual  of  Barometry.  The  method  used  in  obtaining  an  es  for  stations 
not  listed  in  the  table  is  to  select  the  data  from  the  station  in  the 
table  that  has  similar  topographic  and  climatological  conditions  as  the 
station  in  question.  Figure  3.6  shows  the  PAM  II  stations  and  the  NWS 
station  from  which  saturated  vapor  pressure  values  were  obtained  for 
the  PAM  II  stations.  The  hatched  regions  indicate  the  PAM  II  stations 
for  which  data  was  obtained  by  averaging  data  from  two  NWS  stations. 

To  illustrate  the  somewhat  involved  process  of  pressure  reduction 
described  above,  suppose  we  wish  to  obtain  sea-level  pressures  for  1200 
UTC  25  February  1986.  First,  the  temperature  argument  for  each  PAM  II 
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pressure  data  for  the  PAM  II  stations  for  use  in  the  Sea-Level  Pressure 
Reduction  Procedure  outlined  in  the  Manual  of  Barometry  (1961).  The 
unhatched  areas  indicate  that  the  data  for  the  PAM  II  stations  were 
obtained  from  the  lone  NWS  station  in  the  area. 


station  is  calculated  using  the  temperatures  recorded  at  1200  UTC  and 
0000  UTC.  Based  on  the  temperature  argument  for  each  station,  a 
saturated  vapor  pressure  is  found  from  Table  7.2.1  in  the  Manual  of 
Barometry.  Next,  the  standard  lapse  rate  and  humidity  corrections  are 
found  from  Table  7.3  in  the  Manual  of  Barometry  and  F  is  found  using 
Table  7.4.1  in  the  Manual  of  Barometry.  Next,  the  mean  virtual 
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temperature  is  calculated  from  equation  (3.7). 

Once  the  mean  virtual  temperature  is  found,  r  can  be  calculated 
from  equation  (3.6)  with  K  =  0.0266895‘R/m.  Sea  level  pressure  is  then 
found  by  multiplying  the  corrected  station  pressure  by  r. 

The  resulting  sea-level  pressures  for  the  PAM  II  stations  for  our 
example  are  plotted  along  with  the  NWS  sea-level  pressures  and  the 
resulting  analysis  is  shown  in  Figure  3.7a.  This  analyzed  pressure 
field  matches  well  with  the  National  Meteorological  Center  (NMC) 
analysis  of  the  same  time  (Figure  3.7b).  Since  the  analyses  in  Figure 
3.7  agree  then  it  follows  that  the  PAM  II  pressures  agree  well  with  the 
NWS  pressure  field. 

This  method  of  reduction  to  sea  level  is  rather  complicated  as  one 
must  consult  tables  to  obtain  the  parameters  and  F  in  (3.7).  In  an 
effort  to  obtain  a  more  simplified  method,  an  attempt  was  made  to 
determine  the  magnitude  of  the  virtual  temperature  deviation 
permissible  in  order  to  maintain  a  sea-level  pressure  within  0.3  mb. 
Figure  3.8  shows  the  change  in  virtual  temperature  required  to  obtain  a 
sea-level  pressure  deviation  of  0.3  mb  for  station  elevations  up  to 
800  m.  From  Figure  3.8  it  can  be  seen  that  the  virtual  temperature 
deviation  decreases  with  height  to  around  1.5  *F  at  800  m.  The  value 
of  F  in  Equation  (3.7)  is  usually  around  10’F  so  based  on  Figure  3.8,  F 
must  be  retained.  The  value  of  the  term  (aH/2  +  e^)  is  between  1.5'F 
and  2.5*F.  This  term  begins  to  become  significant  for  stations  with 
elevations  near  600  m.  The  entire  term  can  be  approximated  by  a  value 
of  2*F  because,  even  at  the  term's  extremes,  the  virtual  temperature 
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-  t  +  t  +  F  +  2*F. 
TiV  *  T0  +  t, 


(3.9) 


EQliatlon  ,3.9)  still  contains  F  which  retires  the  use  of  Table 

.  „063)  Using  Table  7.4.1  (reproduced 

7.4.1  In  the  Manual  of  Barometry 

in  Figure  3.9)  an  empirical  formula  «•»  derived  to  so  ve 

,0ll0W5:  F,t )  per  -F  for  a  t,  of  0-F  «•  calculated  and  found 

A„  average  F,V  ~  ^  ^  ^  _  cemper.ture 

-  *  °'4976'  ^  *Ven;" ‘7 „  .as  calculated  and  found  to  be 

«-»  0Vet  *  11  r“86  ,  for  F  w&8  developed  using  the  F,V  for  * 

0.434.  An  empirical  formula  for  F  was 


t,  of  0*F  and  a  tn  of  35*F  from  Table  7.4.1  as  the  base  value.  The 
empirical  formula  Is 

F  =  14.3’F  +  0.4375(tn  -  35*F)  -  0.434t,.  (3.10) 

Calculating  a  sea  level  pressure  from  equations  (3.9)  and  (3.10) 
and  comparing  the  results  Illustrated  In  Figure  3.7  reveals  that  these 
new  sea  level  pressures  are  within  0.3  mb  for  those  stations  with 
elevations  of  305  m  or  less.  Therefore,  it  can  be  concluded  that,  to  a 
first  approximation,  this  revised  method  can  be  used  to  reduce  the  PAM 
II  station  pressures  to  sea  level.  Table  3.8  illustrates  the 
comparison  of  the  two  methods  for  1200  UTC  25  February  1986. 
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TABLE  7.4.1 

Correction  for  Plateau  Effect  and  Local  Lapse  Rate  Anomaly,  F(t,), 
Continental  US.  Stations  Having  Elevations  of  305  gpm 
(1,000  Feet)  or  Lower 
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Figure  3.9.  Table  7.4.1  from  the  Manual  of  Barometry. 
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Table  3.8  Comparison  between  using  the  Manual  of  Baroaetry  (MoB) 
method  and  the  revised  method  of  reducing  PAM  II  station  pressures  to 
sea  level  for  1200  UTC  25  February  1986.  Units  of  the  sea  level 
pressures  are  millibars. 
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4.  ALTERNATE  METHOD  FOR  PRESSURE  REDUCTION 

The  NWS  reduction  procedure  is  rather  involved  and  includes  a 
number  of  assumptions  to  arrive  at  the  sea-level  pressure.  In  this 
section,  another  method  of  pressure  reduction,  introduced  by  Pielke  and 
Cram  (1987),  i6  explored  to  eliminate  the  arbitrariness  of  the  standard 
reduction  procedure.  The  method  first  calculates  a  geostrophic  wind 
field  and  then  calculates  the  surface  pressure  corresponding  to  this 
wind  field. 

In  this  research,  both  the  standard  pressure-reduction  method  and 
the  Pielke  and  Cram  method  will  be  applied  to  a  cold-air  damming 
situation  on  25  January  1986  (GALE  IOP  2).  The  1200  UTC  data  are 
being  used  because  in  the  routine  National  Meteorological  Center  (NMC) 
analysis,  the  observed  surface  winds  are  nearly  perpendicular  to  the 
isobars  over  North  and  South  Carolina.  The  pressure  and  geostrophic 
wind  fields  computed  from  each  method  will  be  compared  in  order  to 
determine  the  effect  of  the  assumptions  in  the  standard  reduction 
procedure.  The  geostrophic  wind  fields  calculated  from  each  method 
will  be  compared  to  the  observed  surface  wind  field  in  an  attempt  to 
determine  the  true  ageostrophic  wind  components  during  a  cold-air 
damming  situation. 

4.1  Theory 

The  method  introduced  by  Pielke  and  Cram  (1987)  uses  a  terrain- 
following  coordinate  system  to  obtain  the  horizontal  pressure  field 
from  observations  obtained  on  a  sloping  surface.  The  vertical 
coordinate  of  this  system  was  developed  by  Pielke  and  Martin  (1981)  and 
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is  defined  as: 


Z  -  Zg 

a  =  s  -  .  (4.1) 

S  -  ZG 

where  o  is  the  generalized  vertical  coordinate,  s  is  an  arbitrary 
height  in  the  atmosphere  (usually  taken  as  the  top  of  the  model),  Zg  is 
the  terrain  elevation,  and  z  is  the  height  above  the  surface  (Pielke 
and  Martin,  1981).  Figure  4.1  illustrates  an  example  of  the  terrain 
following  coordinate  system  (Pielke,  1984). 


Figure  4.1.  Schematic  of  the  terrain-following  coordinate  system,  o 
(from  Pielke,  1984). 


Pielke  (1984)  shows  that  the  simplified  momentum  equation  for 
terrain-following  flow  can  be  written  as: 
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where  u  and  v  are  the  o-parallel  flow  in  the  east-west  (x)  and  north- 
south  (y)  directions,  respectively,  f  is  the  Coriolis  parameter,  ©  is 
the  potential  temperature,  g  is  the  acceleration  due  to  gravity, 
n  =  CpT/0  =  c?(p/p0)^cp,  and  T  is  the  temperature.  The  scaled 
pressure,  n,  is  directly  related  to  pressure,  p.  The  pressure-gradient 
term  is  now  written  as  the  sum  of  the  n-gradient  along  a  sigma  surface 
and  the  terrain -gradient  term  (Pielke  et  al.,  1985). 

Neglecting  acceleration,  the  momentum  equation  (4.2)  reduces  to 
the  geostrophic  wind  equation  in  the  o-coordinate  system: 


0  9tc  g  a  -  sdz G 

v„  =  -?■ - - - . 


dx 


dx 


0  dn  i  g  a  -  sdzG 

U  n  /•  '  "t"  /• 


dy 


dy 


(4.3) 


On  a  horizontal,  level  surface  (no  terrain  slope)  Equation  (4.3) 
reduces  to  the  geostrophic  wind  in  the  x-y-z  coordinate  system: 

0  dn 

Vg-T  ’ 

3x  (4.4) 

0  dn 


At  the  surface,  z  =  z^,  so  from  (4.1)  o  =  0  and  (4.3)  becomes: 


0  dn  g  dzc 

v  =T  — +  f - 

‘  dx  ‘  dx 
0  dn  g  dzG 

Ue=_Ti7_T'a7 ' 


(4.6) 


This  equation  for  the  surface  geostrophic  wind  is  then  used  to  obtain 
the  horizontal  pressure  gradient  over  flat  terrain: 
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0  9tc  g  9zG  e  dJt 

Vg  f  ax +  f  dx  ~Ta7  ’ 

__  0  dn  _g  9zg__  0^071  (4,6) 

8  f  9y  f  9y  f  By 

The  pressure  on  a  flat  terrain  surface,  ff,  can  be  obtained  from 
the  equation: 


0  7t  0  TC  0  .vgf  9  /U2^s 

+  =_(_!_) - (  s_) 

0x2  9y“  9x  0  0y  0 


(4.7) 


Here.  u?  and  v?  are  obtained  from  (4.5). 

The  component  of  the  o-system  pressure  gradient  of  interest  is  the 
one  on  a  flat  surface.  Since  in  a  flat,  mean  sea-level,  z-system  the 
geostrophic  wind  is  non-divergent  except  for  the  variation  of  f,  v?  is 
differentiated  with  respect  to  x  and  u?  is  differentiated  with  respect 
to  y  and  subtracted  to  yield  an  elliptic  equation  for  n  (Pielke  and 
Cram.  1987).  Note  that  equation  (4.7)  is  the  geostrophic  relative 
vorticlty  equation  in  this  coordinate  system. 

4.2  PROCEDURE 

The  data  involved  in  this  procedure  are  obtained  from  the  stations 
in  Figure  3.4  and  the  fifty  PAM  II  stations.  First,  the  potential 
temperature,  6.  and  the  modified  pressure  argument,  n,  are  computed 
from  the  temperature  and  station  pressure.  Then,  a  gridded  array  of 


height,  e,  and  n  is  obtained  from  the  Barnes  objective  analysis  with  a 
weight  parameter,  k,  of  10000  km*  and  a  convergence  parameter,  r,  of 
0.3  (same  as  discussed  in  Chapter  3). 

The  three  gridded  arrays  are  used  in  obtaining  the  u  and  v 
components  of  the  geostrophic  wind  for  each  grid  point  from  equation 
(4.5).  From  the  u  and  v  components,  a  geostrophic  wind  field  is 
generated  and  is  shown  in  Figure  4.2. 

The  final  step  is  obtaining  the  pressure  field  that  would  produce 
the  geostrophic  wind  components  calculated  in  the  previous  step.  To 
obtain  the  pressure  field,  equation  (4.7)  is  solved  using  a  relaxation 
method.  The  right  side  of  equation  (4.7)  is  called  the  forcing 
function.  A  first  guess  of  the  ff  field  is  provided  and  the  computer 
corrects  this  guess  until  the  resulting  ff  field  is  within  a  certain 
tolerance  of  the  forcing  function  (i.  e.  the  ff  field  is  relaxed). 

For  this  research,  the  tolerance  used  is  0.01  *K*1s'1  and  the 
relaxation  factor  is  0.30;  however,  later  experiments  indicate  a 
relaxation  factor  around  0.40  is  more  efficient.  For  the  boundary 
conditions,  Pielke  and  Cram  (1987)  employed  pressures  reduced  to  sea 
level  using  a  standard  lapse  rate.  As  a  result,  this  research  employs 
altimeter  settings  for  the  boundary  conditions  and  a  Barnes  analyzed 
field  of  altimeter  settings  for  the  first  guess. 

The  altimeter  settings  are  converted  to  ff  using  the  equation: 

5  =  Cj  (p/1000)!  ,  (4.8) 

with  K  =  R/cf  and  p  replaced  with  the  altimeter  setting.  Once  the  ff 
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Figure  4.3.  Horizontal  pressure  field  corresponding  to  the  geostrophic 
wind  field  shown  In  Figure  4.2. 


4.3  RESULTS 

The  surface  pressure  analysis  of  Figure  4.3  has  the  same  general 
pattern  as  the  sea-level  pressure  analysis  shown  In  Figure  4.4.  A  high 
pressure  ridge  extends  through  the  Carolinas  with  an  inverted  trough 
along  the  coast.  However,  the  ridge  over  western  North  and  South 
Carolina  appears  to  be  shifted  to  the  west  in  Figure  4.3.  The  standard 
pressure  reduction  procedure  resulted  In  the  sea-level  pressures  being 


Figure  4.4.  Sea-level  pressure  analysis  for  1200  UTC  25  January  1986. 

about  1  mb  higher  than  the  surface  pressure  field  calculated  in  the 
previous  section.  The  pressure  gradient  over  the  mountains  appears  to 
be  greater  in  Figure  4.3  than  in  Figure  4.4.  In  a  cold-air  damming 
situation,  one  would  expect  a  tight  packing  of  the  isobars  at  the 
mountains  due  to  the  cold  air  being  dammed  against  the  mountains.  It 
would  seem  that  the  assumptions  in  the  standard  reduction  procedure 
reduce  this  damming  effect  on  the  pressure  field  and  result  in  a  slight 
underestimate  of  the  gradient. 
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One  other  item  to  note  about  the  comparison  of  Figures  4.3  and  4.4 
is  the  isobars  over  the  coastal  waters  do  not  match  as  would  be 
expected.  In  both  procedures,  the  pressure  data  recorded  at  the  marine 
stations  did  not  need  to  be  reduced,  but  the  altimeter  equation 
subtracts  0.3  mb  from  the  station  pressure.  This  subtraction 
accounts  for  the  difference  in  the  isobars  over  the  water. 

Moving  on  to  the  comparison  between  the  surface  geostrophic  wind 
field  depicted  in  Figure  4.2  and  the  geostrophic  wind  field  calculated 
from  the  sea-level  pressures  depicted  in  Figure  4.5,  the  same  general 
pattern  is  found;  however,  the  ridge  in  the  wind  field  is  shifted 
westward  in  Figure  4.2.  The  cause  of  this  shift  could  be  the  result  of 
the  dominance  of  the  terrain  term  in  equation  (4.5).  In  the  absence  of 
a  pressure  gradient,  the  effect  of  the  terrain  is  to  produce  a  ridge 
over  the  mountains  and,  on  a  flat  terrain,  the  effect  of  the  pressure 
term  is  to  produce  a  trough  over  the  mountains.  These  two  effects 
should  almost  cancel  each  other  in  equation  (4.5),  but  because  high 
pressure  is  dominating  the  region,  the  contribution  of  the  pressure 
term  is  reduced  allowing  the  westward  shift  of  the  ridge  in  the  wind 
field  in  Figure  4.2. 

Overlaying  Figure  4.2  over  Figure  4.3  reveals  that  there  is  cross- 
isobaric  flow  in  the  north-central  and  northwest  portions  of  the  map. 
According  to  the  geostrophic  wind  balance,  since  Figure  4.3  is 
determined  by  the  geostrophic  wind  field  (Figure  4.2)  the  winds  should 
be  parallel  to  the  isobars.  A  number  of  sample  calculations  were 
accomplished  to  determine  any  errors  in  the  computer  result.  No 
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Figure  4.5.  Geostrophic  wind  field  calculated  using  the  sea-level 
pressures  for  1200  UTC  25  January  1986. 

problems  were  found,  so  it  would  appear  that,  in  the  relaxation 
process,  the  forcing  function  defined  by  the  geostrophic  wind  field  is 
attempting  to  force  the  pressure  field  to  conform  as  is  expected,  but 
the  boundary  conditions  Imposed  on  the  pressure  field  are  not  allowing 
the  pressure  field  to  conform  completely.  Indeed,  in  comparing  the  two 
pressure  fields  (Figures  4.3  and  4.4)  there  is  a  slight  westward  shift 
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in  the  ridge  axis  in  the  surface  pressure  field  suggesting  the  surface 
pressure  field  is  attempting  to  conform  to  the  geostrophic  wind 
balance.  A  further  complication  in  this  case  is  that  the  mountains  are 
located  along  the  boundary  adding  the  effect  of  the  mountains  to  the 
effect  of  the  boundary  conditions.  In  the  Pielke  and  Cram  (1987)  study 
a  much  larger  area  was  used  and  there  was  evidence  of  cross-isobaric 
flow  towards  the  boundary  of  the  area.  Therefore,  had  a  larger  area 
been  employed  in  this  research,  the  pressure  field  could  be  expected  to 
conform  more  totally  to  the  geostrophic  wind  field. 

The  final  item  to  address  in  this  section  is  whether  the  method 
for  calculating  geostrophic  winds  described  in  this  section  produces 
less  of  an  ageostrophic  component  when  compared  to  the  observed  wind 
field  than  the  geostrophic  wind  calculated  with  the  sea-level 
pressures.  Figure  4.6  shows  the  wind  vectors  of  the  observed  winds 
(o),  the  surface  geostrophic  winds  from  Figure  4.2  (SG),  and  the  sea- 
level  geostrophic  winds  from  Figure  4.5  (G)  for  selected  points.  From 
Figure  4.6,  a  significant  improvement  can  be  seen  in  SG  over  G  in  the 
western  Piedmont.  The  surface  geostrophic  winds  are  much  closer  to  the 
observed  winds  than  the  geostrophic  as  was  hoped.  However,  over  the 
rest  of  the  area,  there  is  a  slight  improvement  in  direction  of  SG  over 
G,  but  not  much  improvement  in  the  speed.  The  observed  winds  are  still 
ageostrophic,  but  SG  seems  to  be  a  better  representation  of  the 


observed  winds  than  G. 


Figure  4.6.  Plot  of  wind  vectors  of  the  observed  winds  (0),  surface 
geostrophic  winds  (SG),  and  the  sea-level  geostrophic  winds  (G)  for 
selected  points  at  1200  UTC  on  25  January  1986. 


5.  CONCLUSIONS 


The  basic  objective  of  this  research  is  to  enable  the  pressure 
data  recorded  at  the  PAM  II  stations  during  GALE  to  be  used  in 
mesoscale  analyses.  This  objective  was  to  be  met  by  first  correcting 
the  pressures  and  then  reducing  the  corrected  pressures  to  sea  level. 

Pressure  corrections  were  found  for  each  of  the  PAM  II  stations 
using  the  Barnes  objective  analysis  scheme  to  interpolate  an  altimeter 
setting  to  the  PAM  II  stations  from  the  altimeter  settings  reported  by 
the  NWS  stations.  The  Barnes  analysis  scheme  was  selected  over  direct 
one-to-one  comparisons  between  the  PAM  II  stations  and  their  nearest 
NWS  neighbor  for  two  reasons.  First,  not  every  PAM  II  station  was 
located  near  a  NWS  station.  Second,  the  Barnes  analysis  was  expected 
to  be  a  better  method  because  it  incorporated  data  from  several 
stations  into  the  interpolated  altimeter  setting.  However,  for  the 
thirteen  cases  Included  in  the  comparisons,  a  Chi-squared  test  showed 
that  the  Barnes  analysis  did  not  provide  a  statistically  significant 
improvement  over  the  direct  one-to-one  comparison,  or  buddy  check. 

The  interpolated  altimeter  setting  was  compared  to  the  altimeter 
setting  calculated  from  the  observed  station  pressure  for  each  PAM  II 
station  to  obtain  an  offset.  The  offsets  from  thirteen  cases  were 
averaged  to  obtain  the  pressure  corrections  for  each  station.  The 
corrections  were  applied  to  the  station  pressure  and  the  corrected 
pressures  were  reduced  to  sea  level  using  the  method  described  in  the 
Manual  of  Barometry  (1963).  The  analysis  of  the  reduced  PAM  II 
pressures  and  the  reported  NWS  pressures  matched  well  with  the  NMC 
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analysis  for  two  separate  cases,  suggesting  that  the  pressure 
corrections  are  valid  for  the  GALE  period. 

The  procedure  outlined  in  the  Manual  of  Barometry  (1963)  does  not 
lend  itself  easily  to  manipulation  by  a  computer  as  two  of  the 
corrections  to  the  mean  virtual  temperature  involve  the  use  of  tables. 

The  result  is  a  procedure  that  takes  a  great  deal  of  time  to  find  sea- 
level  pressure.  An  empirical  formula  for  calculating  the  mean  virtual 
temperature  without  entering  tables  was  developed  to  allow  a  computer 
to  calculate  the  sea-level  pressure  for  each  station.  The  sea-level 
pressures  calculated  from  this  empirical  formula  were  within  0.3  mb  of 
the  pressures  computed  using  the  tables.  So,  to  a  first  approximation, 
this  formula  can  be  used  in  calculating  sea-level  pressures. 

While  the  first  two  objectives  have  been  wet,  the  sea-level 
reduction  procedure  still  involves  a  number  of  assumptions  in  the  lapse 
rate  through  a  fictitious  column  of  air  which  was  shown  in  previous 
work  to  create  errors  in  the  sea-level  pressure  analysis.  A  final 
objective  of  this  research  is  to  study  an  alternate  method  of  pressure 
reduction.  This  method  produces  a  horizontal  pressure  field  from 
observations  on  a  sloping  surface  by  first  calculating  the  geostrophic 
wind  components  defined  on  a  terrain-following  coordinate  system.  In 
this  system,  the  pressure  gradient  term  in  the  geostrophic  wind 
equation  consists  of  a  modified  pressure  gradient  term  plus  a  height 
gradient,  or  slope,  term.  From  the  geostrophic  wind  components,  the 
corresponding  horizontal  pressure  field  is  obtained  through  the  use  of 


a  relaxation  method. 
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The  resulting  surface  pressure  field  (Figure  4.3)  has  the  same 
general  pattern  as  the  sea-level  pressure  analysis  (Figure  4.4),  but 
there  appears  to  be  a  stronger  gradient  in  the  surface  pressure  field 
over  the  mountains  than  the  sea-level  pressure  analysis.  During  cold- 
air  damming  situations,  a  tighter  packing  of  the  isobars  against  the 
mountains  would  be  expected  as  colder  air  is  being  dammed  up  against 
the  mountains.  The  sea-level  pressure  values  are  about  1  to  2  mb 
stronger  than  those  obtained  from  the  surface  pressure  analysis.  It 
would  appear  that  the  assumptions  used  in  the  standard  sea-level 
pressure  reduction  method  reduce  the  effect  off  the  cold-air  damming 
resulting  in  a  weaker  gradient  over  the  mountains  and  lead  to 
artificially  higher  pressure. 

The  surface  geostrophic  winds  are  not  parallel  to  the  surface 
pressure  isobars  as  is  expected.  In  the  relaxation  process,  the 
forcing  function  defined  by  the  geostrophic  wind  field  attempts  to 
create  a  pressure  field  so  the  winds  are  parallel  to  the  isobars,  but 
the  boundary  conditions  imposed  on  the  pressure  field  along  with  the 
mountains  located  at  the  western  boundary  may  prevent  the  pressures 
from  conforming  completely.  Using  a  larger  area  so  the  mountains  are 
away  from  the  boundary  should  produce  the  expected  results. 

The  surface  geostrophic  winds  calculated  from  equation  (4.5)  were 
expected  to  be  an  improvement  over  the  geostrophic  winds  calculated 
using  the  sea-level  pressures.  The  observed  winds  were  expected  to  be 
less  ageostrophlc  relative  to  the  surface  geostrophic  winds  than  to  the 
sea-level  geostrophic  winds.  Figure  4.6  shows  there  to  be  a  direction 
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improvement,  but  not  much  of  a  speed  improvement  over  much  of  the  area 
of  interest;  however,  over  the  western  Piedmont,  where  the  winds  are 
most  ageostrophic,  the  surface  geostrophic  winds  are  much  closer  to  the 
observed  winds  than  the  geostrophic  winds.  One  reason  there  is  not 
much  improvement  in  the  speed  elsewhere  might  be  that  the  observed 
winds  were  so  light  (around  5  m/s).  Even  so,  it  would  appear  that  the 
method  Introduced  by  Pielke  and  Cram  (1987)  is  a  better  representation 
of  the  surface  geostrophic  winds  than  using  the  sea-level  pressures. 
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6.  FUTURE  RESEARCH 

Further  research  is  needed  to  continue  to  validate  the  pressure 
corrections  derived  in  this  research.  This  could  be  done  by  applying 
the  pressure  corrections  to  other  cases  during  GALE.  One  such  case 
would  be  the  cyclogenesis  that  occurred  on  13  March  1986.  The  pressure 
corrections  could  be  applied  to  the  5-minute  data  to  determine  if  the 
corrected  pressure  data  permit  the  detection  of  an  important,  but  weak 
and  short-lived  low-pressure  center  that  developed  in  the  PAM  II 
region. 

Other  applications  of  these  pressure  corrections  would  be  the 
study  of  mesoscale  thunderstorm  cases,  such  as  the  one  depicted  in 
Figures  3.2  and  3.3  to  determine  whether  a  trigger  can  be  detected  in 
the  mesoscale  pressure  pattern.  Finding  a  trigger  to  these  storms 
would  enable  better  prediction. 

Future  research  could  also  be  directed  toward  studying  the 
empirical  formula  derived  in  Chapter  3  for  reducing  station  pressure  to 
sea-level  in  order  to  create  a  formula  for  outside  the  GALE  region  and 
GALE  period. 

More  important  than  working  on  the  empirical  formula  would  be 
continued  testing  of  the  alternate  method  of  pressure  reduction 
introduced  by  Pielke  and  Cram  (1987)  in  order  to  determine  whether  the 
method  can  be  used  outside  the  Great  Plains.  One  such  test  could  be  on 
a  stronger  cold-air  damming  case  than  was  presented  in  this  research  to 
determine  if  the  method  produces  winds  closer  to  the  observed  winds 
over  a  larger  area. 


In  computing  the  n  field  from  the  geostrophic  wind,  a  much  larger 
area  needs  to  be  employed  to  ensure  the  terrain  feature  is  away  from 
the  boundary  and  the  boundary  conditions  themselves  do  not  interfere 
with  the  computation  of  the  2  field  in  the  area  of  interest;  in  this 
case  the  GALE  region. 
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APPENDIX  8.1 

EARLY  METHODS  OF  SEA-LEVEL  PRESSURE  REDUCTION 

The  first  method  employed  by  the  Weather  Service  was  used  from 
1870  to  mid- 1881  and  employed  Guyot's  Table  XIX  (Bigelow,  1900).  This 
table  provided  the  height,  in  feet,  of  a  column  of  air  corresponding  to 
0.10  in  of  mercury  at  different  temperatures  and  elevations.  The 
procedure  was  to  increase  the  temperature  of  the  higher  elevation 
station  to  that  of  a  lower  station  by  using  a  standard  lapse  rate. 
Entering  the  table,  one  then  could  obtain  the  number  of  feet  per  0.10 
in  mercury  for  each  station.  Next,  one  found  the  mean  of  these  two 
heights  and  divided  the  elevation  difference  by  this  mean  to  obtain  the 
correction  to  be  added  to  the  pressure  of  the  higher  elevation  station 
(Bigelow,  1900). 

This  procedure  was  found  to  have  involved  two  procedures  which 
greatly  exaggerated  the  diurnal  variation  of  the  isobars  making  them 
almost  useless  for  forecasting  (Bigelow,  1900).  The  first  problem  was 
that  the  temperature  argument  was  taken  as  the  reported  surface 
temperature  at  the  hour  of  observation  rather  than  for  a  time  average 
for  the  plateau  column.  Thus,  the  whole  air  column  was  considered 
subject  to  a  diurnal  range  of  temperatures  which  in  fact  is  only 
characteristic  of  the  air  near  the  surface.  The  second  problem  was 
that  Table  XIX  was  only  entered  once  to  get  the  mean  reduction  number 
instead  of  twice;  once  to  reduce  the  higher  station  to  the  elevation  of 
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the  lower  station  and  once  to  reduce  the  lower  station  to  sea  level 
(Bigelow,  1900). 

In  an  effort  to  eliminate  the  problems  encountered  using  Guyot’s 
tables,  the  Weather  Service  adopted  a  method  that  employed  monthly 
constants  for  the  reduction  of  pressure  developed  by  Professors  Abbe 
and  Upton  (1882).  The  Abbe-Upton  method  took  three  steps  in 
determining  the  observed  station  temperature.  The  first  was  to  compute 
the  latitudinal  average  temperature  at  sea  level.  The  second  step  was 
to  calculate  an  average  lapse  rate  for  a  given  latitude.  The  final 
step  was  to  calculate  the  variation  of  lapse  rate  along  a  given 
latitude  (Abbe  and  Upton,  1882).  The  values  of  the  temperature  and 
lapse  rate  were  used  to  calculate  a  mean  monthly  temperature  of  the 
station  (Abbe  and  Upton,  1882).  This  mean  monthly  temperature  was  then 
used  in  Upton's  formula  for  pressure  reduction. 

Three  problems  were  soon  found  with  this  method.  The  first  was 
that  the  sea  level  reduction  of  pressure  ought  to  involve  two  other 
corrections,  one  depending  on  the  deviation  of  pressure  from  the  mean 
and  the  other  depending  on  the  deviation  of  temperature  from  the  mean. 
These  deviations  occur  almost  constantly  and  must  be  taken 
into  consideration,  but  are  not  considered  when  using  monthly 
constants.  The  second  problem  was  that  Upton's  computed  lapse  rates 
were  greater  In  the  winter  than  summer,  whereas  actual  observations 
showed  them  to  be  greater  in  the  summer.  The  final  problem  is  that  the 
method  did  not  produce  useful  pressure  fields  (Bigelow,  1900). 

As  a  result  of  the  failure  of  the  monthly  constants,  William 
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Ferrel  was  commissioned  by  the  Weather  Service  to  develop  a  new  method 
for  the  reduction  of  pressures.  Ferrel  (1886)  made  four  improvements 
to  the  earlier  methods.  The  first  was  the  substitution  in  the 
reduction  procedure  of  the  mean  diurnal  temperature  instead  of  the 
temperature  at  the  hour  of  observation.  The  second  improvement  was  the 
introduction  of  a  mean  lapse  rate  to  increase  the  surface  temperature 
as  an  argument  to  the  mean  temperature  of  the  air  column.  Ferrel 
(1886)  adopted  a  lapse  rate  of  1*F  per  600  feet  and  applied  it  for  half 
the  elevation  of  the  station.  The  third  improvement  was  the 
introduction  of  a  thermal  pressure  term  which  Ferrel  defined  as  CAH; 
where  C  is  a  constant  with  a  value  of  0.00105  in/’Fft,  A  is  the 
departure  of  the  temperature  from  the  normal,  and  H  is  the  station 
elevation  in  feet  (Ferrel,  1886).  The  last  improvement  was  the 
introduction  of  a  plateau  correction.  Ferrel  found  that  all  the 
reduced  pressures  for  all  high  stations  were  too  great  in  January  as  a 
result  of  using  temperatures  that  were  too  cold,  and  too  low  in  July  as 
a  result  of  using  temperatures  that  were  too  warm.  It  was  necessary  to 
apply  some  correction  to  the  reductions  for  these  extremes.  The 
correction  was  determined  on  the  principle  that  the  differences  of  the 
reduced  pressures  of  the  high  stations  in  January  and  July  should  not 
be  greater  than  those  of  the  low-level  stations,  which  Ferrel  found  to 
be  0.15  in.  (Ferrel,  1886).  The  plateau  correction  Ferrel  developed  is 
defined  as  follows: 


a  =  c  +  v  \H 


(1.1) 
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with  a  =  correction  of  pressure  reduced  to  sea  level  (in). 

c  =  average  departure  of  pressure  at  sea  level  from  the  normal 
annual  pressure  (in). 

CAH  =  is  the  thermal  pressure  term  as  defined  above. 

Ferrel  assumed  c  to  be  a  constant  and  equal  to  0.073  In.  However, 
Bigelow  (1900)  reported  that  c  is  not  a  constant,  but  has  an  annual 
range  from  +0.073  to  -0.073.  The  value  should  be  calculated  monthly  as 
well  as  for  the  mean  for  January  and  July  (Bigelow  (1900). 

In  his  report  of  1900,  Bigelow  introduced  improvements  to  Ferrel's 
techniques  by  adopting  a  variable  lapse  rate  for  the  fictitious  air 
column  and  a  humidity  correction  as  a  function  of  the  surface 
temperature,  and  by  redefining  the  temperature  argument. 

Bigelow  believed  that  the  assumed  lapse  rate  in  the  fictitious  air 
column  should  vary  with  the  location  of  the  station,  the  height  above 
sea  level,  and  the  month.  If  the  assumed  lapse  rate  Is  allowed  to  vary 
in  an  appropriate  manner  with  the  station  temperature  argument 
depending  on  the  station,  it  would  not  be  necessary  to  regard  the 
monthly  temperature  as  anomalous  (Manual  of  Barometry,  1963). 

In  his  method  for  dealing  with  the  humidity  term,  he  took  the 
vapor  pressure  at  the  surface  and  employed  the  following  formula  to 
compute  an  average  vapor  pressure  of  the  column: 

e  =  e/*17  (1.2) 

where  h  Is  the  height  in  meters,  e0  is  the  vapor  pressure  at  sea  level, 
and  e  is  the  vapor  pressure  at  height,  h  (Bigelow,  1900).  Bigelow 
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added  a  correction  to  the  above  formula  for  abnormal  conditions.  This 
formula  is: 

... ,■***'*<  (1.3) 

Finally,  Bigelow  defined  the  station  temperature  argument  as  the 
average  of  the  current  observed  temperature  and  the  observed 
temperature  twelve  hours  previously  (Manual  of  Barometry,  1963). 

The  Weather  Service  adopted  a  method  that  Includes  the  work  of 
Ferrel  and  Bigelow.  The  pressure  reduction  is  accomplished  using  the 
hypsometric  equation  with  the  temperature  replaced  by  the  mean  virtual 
temperature  of  the  air  column.  The  calculation  of  the  mean  virtual 
temperature  includes  the  temperature  argument  defined  by  Bigelow,  a 
standard  lapse  rate  correction,  a  humidity  correction  first  introduced 
by  Bigelow,  and  a  plateau  correction  developed  by  Ferrel  (Manual  of 
Barometry,  1963). 
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APPENDIX  8.2 


This  table  presents  the  PAN  II  station  elevations  used  in  this 
research.  Elevations  are  in  aeters. 


Station 

Elevation 

Station 

Elevation 

01 

2 

26 

271 

02 

18 

27 

373 

03 

50 

28 

722 

04 

93 

29 

764 
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APPENDIX  8.3  (continued) 
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